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We simulated the growth of 2D islands with 2 kinds of diffusion atoms using the kinetic Monte- 
Carlo (kMC) method. As a result, we found that the slow atoms tend to create nuclei and determine 
the island volume distribution, along with additional properties such as island density. We also 
conducted a theoretical analysis using the rate equation of the point-island model to confirm these 
results. 
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Nucleation is an important step in several phenom¬ 
ena, such as growth of clouds, polymers, and crys¬ 
tals. In many cases, heterogeneous nucleation occurs 
on surfaces and at interfaces, and it is also facili¬ 
tated by the presence of impurities. The heterogeneous 
two-dimensional (2D) nucleation is particularly impor¬ 
tant because of its relevance to the epitaxial growth of 
single crystal semiconductors, of which quantum dots 
(QDs)[ 3-Q are one of the most attractive classes. 

Many theoretical and experimental works have been 
conducted on the 2D nucleation and island size- 
distribution(B-[l^ However, most of the studies only 
dealt with one kind of diffusion atoms. We have investi¬ 
gated the volume distribution of QDs grown by molecu¬ 
lar beam epitaxy (MBE), and we are now interested in 
ternary alloy semiconductor QDs, especially in InAlAs 
QDs, which have a luminescence wavelength of 750 nm, 
at which the Si photodetectors exhibit the maximum 
sensitivity (2,1^. InAlAs QDs are ternary QDs, to which 
two group in adatoms. In and Al, with extremely dif¬ 
ferent surface diffusion rates, contribute. We previously 
found that investigating the submonolayer 2D island 
growth provides helpful knowledge on the QD growth 
regime[J,y, and the 2D island growth by two kinds of 
largely different adatoms is of academic interest. In this 
paper, we examine the 2D nucleation and island size- 
distribution using both computer simulations and the¬ 
oretical analysis of the point-island model. We found, 
for the first time, that the adatom with the slower dif¬ 
fusion tends to form nucleation centers and determines 
the size distribution of the islands. 


adatom is chosen to diffuse to its nearest-neighbor site 
with a specific diffusion rate, D. At nucleation, the im¬ 
portant parameter is the critical island size, i, such that 
i + 1 aggregate atoms form a nucleus, which is a min¬ 
imum island. Here we used i = 1, which means that, 
if an adatom encounters another adatom as a neighbor, 
they are frozen at their sites and form a new island, re¬ 
sulting in nucleation. These processes are repeated up 
to a fixed coverage, 9, which is the ratio of the number 
of all deposited atoms against the lattice area; typically, 
we set 9 — 0.3. Here, R = D/F characterizes the sim¬ 
ulation for a fixed 9. We then averaged the data from 
approximately 100 simulation runs. 

Regarding InAlAs QD growth, Al has a much smaller 
diffusion coefficient than In, because of the stronger Al- 
As bond in comparison with the In-As bond till. There¬ 
fore, we included two adatoms with different diffusion 
rates, Df and corresponding to fast In and slow Al, 
respectively. We varied Rf = Df/F and Rg = Dg/F 
and observed their effects on the island size distribu¬ 
tion. The flux ratio of these two adatoms represents 
the alloy composition of the dots, for which we used a 
value corresponding to the In : Al = 0.7 : 0.3 ratio used 
in our previous experiment Q. 

In a submonolayer epitaxial growth process such as 
that of homoepitaxial Fe/Fe(100) @ or heteroepitax- 
ial Fe/Cu(100) [l^. and also the InAs/GaAs QD used 
in our study |^, the volume distribution of islands is 
known to show the scaling property 


Ng = 9S- 


‘f’ (I) ■ 
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First, we modeled 2D island growth using the kinetic 
Monte Carlo Method 0, [1| ■ In this model, we simulate 
the motion of atoms on a 400 x 400 square lattice with 
a periodic boundary condition. At each stage, a new 
adatom is deposited on a randomly chosen site with 
a given deposition rate, F, or a previously-deposited 


where Ng is the density of an island of size s. S' is the 
average size, and fi{s/S) is the scaling function, which 
varies with the critical island size, i. 

Fig. 1(a) shows the scaling plot at Rf = 10® with 
varying Rg. Each trendline has only a single peak, 
which translates leftward with decreasing Rg. On the 
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right side of the plot, we can see a tail larger than the 
emprial curve for a single atomQ, illustrated by the 
solid curve. This means that these island distributions 
experience size dispersion in the case of large islands, 
where s/S > 2. Here, we note that this feature is simi¬ 
lar to our previous experiments on InAlAs QDs Q. 

Fig. 1(b) shows the scaling plot at Rg = 10® with 
varying R f . All curves are the same and this insensitiv¬ 
ity to changes in i?/ is also found for other Rg values. 
The scaling function is also equal to the single atom 
case with single R = i?s,illustrated by the Rf = 10® re¬ 
sult. These results imply that Rg determines the scaling 
function regardless of the value oi Rf. 

For the scaling arguments the following is assumed 

N ^ R-^0^-\ (2) 

Ni - R-^e-’' , (3) 


where N = X]s >2 is the total island density, A^i 
is the monomer density, and v are critical 

exponents. 

Fig. 2 shows the island density for cases with (a) 
varying Rg and (b) varying Rf. As shown in Fig. 2(a), 
the dependence of the density on Rg is similar to the 
behavior of a single atom, and is expressed as A^ ~ 
R-Xi OY N Ki R~^' for hxed 9. According to 12|, Xi = 
i/{i 2) and xi = 1/3 (for i = 1), and our results give 
Xi ~ 0.33 for the single atom case, and xi ~ 0.31(i?/ = 
10®), 0.30(i?/ = 10®) for cases featuring two atoms. On 
the other hand, as shown in Fig. 2(b), the density is 
almost independent of i?/ for i?/ 3> Therefore we 
can say that the slow adatom plays almost the same 
role about nucleation as the only diffusion atom in the 
single atom case. 

Fig. 3(a) and Fig. 3(b) show the ’’unsealed” island 
size distribution with fixed Rf and fixed i?s, respec¬ 
tively. As shown in Fig. 3(a), the distribution strongly 
depends on Rg. This behavior is similar to that of a 
single atom with R. On the other hand, as shown in 
Fig. 3(b), the distributions are almost the same. This 



FIG. 1. The island size scaling function with (a) varied Rs 
and fixed Rf, and (b) varied Rf and fixed Rs. They are 
compared with scaling functions for a single atom i = 1 
(solid line) and i = 0 (dashed line). 



FIG. 2. The island density for (a) fixed Rs and varying Rf, 
and (b) fixed Rf and varying Rg. The dashed lines show 
the case of single atom. 


result inplies previous two features with fixed Rg, same 
scaling functions in Fig. 1(b) and same island densities 
in Fig.2(b). The slow adatom determines not only the 
scaling function but also the unsealed island size distri¬ 
bution itself. 

Fig. 4(a) shows the origins of nucleation, i.e., the 
nucleus densities of islands composed of pairs of fast- 
fast, fast-slow, and slow-slow atoms. At the early stage, 
or for the cases of small coverage shown in the inset of 
Fig. 4(a), the number of nuclei containing slow atoms 
does not exceed the number of fast-fast nuclei. However, 
in the late stage, there are much more slow-slow nuclei 
than fast-fast nuclei, and even more fast-slow nuclei. 
At the same diffusion rate limit, i.e., Rf = Rg, the 
number of nuclei containing slow adatoms, both slow- 
slow and fast-slow combinations, must be 1 — (0.7)®, or 
51% of all nuclei in our 0.3 : 0.7 composition. Here 
however, this ratio is over 90% during the late stage. 
This demonstrates that slow monomers tend to be used 
as nuclei. This feature is clearly seen in Fig.4(b). 

Fig. 4(b) shows monomer density for each kind of 
adatom. Except for very early stage, the slow monomer 
density is much larger than that of the fast monomers. 
Even in later stage in this simulation, when there are 
few or no monomers, the density of slow monomers is 
about 10“® atoms/site, and that of fast monomers is 




FIG. 3. The ’’unsealed” island size distribution with (a) 
varied Rs and fixed Rf, and (b) varied Rf and fixed Rs. 
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FIG. 4. (a) Nuclei density for each atom pair, (b) monomer 
density for each atom composition, with Rf = 10®, Rs = 
10®. The insets magnify the early stage. 


about 10“^ ~ 10“® atoms/site on time average. These 
show that nucleation by slow adatoms originates from 
the large number of slow adatoms remaining on the sur¬ 
face. 

To confirm the above finding, we also studied the ki¬ 
netics of island growth processes theoretically, using the 
rate equation approach Q . The rate equation describes 
the evolution of island size using the mean-field the¬ 
ory, which we extended to growth featuring 2 types of 
adatoms. Then, the fast monomer, N(, slow monomer, 
Nf, and island with s atoms, Ng, densities can be writ¬ 
ten as 


^ = Ff- 2K( (iV/)2 - k(n{ 

- K^N^nI - N{ ^ KfNg , (4a) 

= Fg- K^NINI - 2K{{Nlf 

- k(N (, (4b) 

s^2 

^ = N{{Ki,Ng_^ - KfNg) 

+ N^iK:_,Ng_^-K:Ng) , (4c) 


where Fa is the deposition rate of adatom a (a = f, s, 
and Ff + Fg = F), and Kf is the rate of adatom a 
attachment to islands of size s. We assume that ~ 
DaS^, Da is the diffusion rate of adatom a, and the 
exponent, p, differs by the shape or dimension of the 
islands. Dividing the equations by F, they can then be 
rewritten as 

d ^ £ £ 

- RgNgNl - RfN( Y s^Ng , (5a) 

s^2 

-^ = fs- RsNIN{ - 2Rg{Nlf 

- RfN{Nl - RgNl Y > (5b) 

s=2 


-^ = [RfN{ + RgNl] 

X [(s - l)PNg_, - sPNg] , (5c) 

where Ra = Da/F, 9 = Ft is the coverage, and fa = 
Fa/F is the deposition ratio of adatom a. Rewriting in 
terms of the scaled variables: 9 = 
and Nf‘ = Rj fVf, the equations become 


dN( 

d9 


= ff- 2{N{f - N(Nl 
- rN^Nl - N( Y , 


(6a) 


dNf 

d9 


s\2 


= fg- rN^Nl - 2r{Nl) 

- N( - rN^ Y ’ 


(6b) 


dNg 

d9 


= [N{ + riVi«][(s - l)^’fV,_i - sPNg] , (6c) 


where r = Rs/Rf 1- 

We further assumed the point-island model @ , where 
each island occupies only one site. Then, Ka ^ Da is 
independent of s (p = 0), and X(s =2 = 

N. Hence, 
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dNj 

d9 


dNl 

d9 


= ff-2{N{f-N{N^ 
- rNlN( -N(N , 


s\2 


= fg- rN^N( - 2r{Nl) 
- N{ Nl - rN(N , 


dN dNs riCrf r f 




(7a) 

(7b) 

(7c) 


In the early stages (low coverage, 0 <C 1), the number 
of monomers is quite low, and which is larger than that 
of islands, so that N <C fV“. Therefore, the last 4 terms 
in Eqs. 7(a) and 7(b) can be neglected, and therefore 

n{ = ff9 (8a) 

Nl = fg9 (8b) 

N r^{ff + rfg)Rf9^ ffRf9^ . (8c) 

In later stages {9 1), the island density has in¬ 

creased while the monomer density has decreased due 
to shortened diffusion length. Hence, N ^ 7V“ and 
dNi /d9 <C fa- Therefore, 

dN( 


d9 

dNf 

d9 

dN 

d9 


ff-N(N:^0, 

fg - rN^N ~ 0 , 
[N( -briVf][7V/ 


(9a) 

(9b) 


N(] 


(9c) 
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Using the condition r <C 1, // + fs/r = fs/r, the solu¬ 
tion can be rewritten as 

n( , ( 10 a) 

, ( 10 b) 

N^flRjhi. ( 10 c) 


Apart from the factor oi fa, N is equal to the sin¬ 
gle atom case solution (replacing R with Rg), N ^ 
ii-l/30l/3[6], in accordance with Fig. 2. In addition, 
the property n( ^ rN^ <C Nf corresponds to Fig. 
4(b). Moreover, the fact that the island distribution 
Ng is independent of Rf, seen in Figs. 1(b) and 3(b), 
can easily be obtained by using eqs. (5c) and (10a) 
-( 10 c) with p = 0 iteratively. 

Also, we can discuss the contents of the nuclei 
7 Y/-/,]\r/-s, and as shown in Fig. 4 (a). Here 

j\j-a-/3 j- 0 pj- 0 g 0 jits the density of islands which have nu¬ 
cleated using atoms a and j5 {a, fi = f or s, arbitrary 
order). The rate equations of these densities can be 
written as 


dNf-f 

dt 

dNf-^ 

dt 

dt 


K(iN(r , 

k(n(+ K^NIN( , 

KliNff . 


( ll a) 

( ll b) 

( ll c) 


These can be rewritten, using the point-island model, 
as 


dNf-f 

d§ 

dNf-^ 

dO 

dN^-‘‘ 

dO 


= {Nlf 


= (1 -f r)NlNl 


= r(Nlf . 


n( Nl 


( 12 a) 

( 12 b) 

( 12 c) 


In the early stage, from eqs.( 8 a)-( 8 c), the solutions 
are 


= ffRfO^ = ffN , 

(13a) 

= fffgRfd^ = fgN , 

(13b) 

TV®-® = rf^RfO^ = rf!fJ^N . 

(13c) 


These results show that nucleus composition depends 
only on the flux ratio except slow-slow nuclei, which 
have much lower density than others, in the early stage, 
and this is consistent with the results displayed in the 
inset of Fig. 4(a). 

In the late stage, using Eqs.(10) and (12), the solu¬ 
tions are 

^ RJ^Rlfjjrhi = rfjf-^N , (14a) 


Nf-^ ^Rg-^fffje"^ = ffN , (14b) 

~ Rj^flei = fgN . (14c) 


This means that the ratio of the nucleus is ^ : 
Nf~‘^ : and therefore 

N^-f A^-f“®,A^®“® when the difference between // 
and fg is not substantial. This result shows that slow 
adatoms tend to be used as nuclei more frequently than 
fast adatoms, and this agrees with the trend shown in 
Fig. 4(a). The number of s-s nuclei in simulation is not 
so little as eqs.(14), because it includes nuclei in the 
early stage like as eqs.( 12 ). 

In summary, we have simulated submonolayer epi¬ 
taxial growth with 2 kinds of adatoms and reached the 
conclusion that slower adatoms tend to form the nuclei 
of islands and also determine island density and size 
distribution. We also conducted a theoretical analysis 
using the rate equation of the point-island model and 
confirmed these results. 

We believe our finding, that the slow adatom deter¬ 
mines the nucleation and volume distribution, is univer¬ 
sal, and therefore can also be applied to the growth of 
InAlAs quantum dots.It is worth noting a report that 
the results for 3D islands show similar properties to 
those of 2D islands 0 . 
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